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ABSTRACT 

This  report  describes  the  current  re-entry  vehicle 
parametric  study  program,  based  on  the  Allen  and  Eggers 
straight  line  trajectory  assumptions.  The  heating  analysis 
considers  a  nonreacting  cold  wall,  with  the  inclusion  of  a 
moving  boundary  conduction  solution  and  boundary  layer 
transition.  The  quantities  determined  are  ablation  thick¬ 
nesses,  insulation  thicknesses,  vehicle  weight  and  aero¬ 
dynamic  coefficients,  for  varying  shape,  material  and 
trajectory  parameters.  The  computer  output  consists  only 
of  final  integrated  values  rather  than  step-by-step  computa¬ 
tions  along  the  trajectory. 
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DEFINITION  OF  SYMBOLS 


Symbol 

A 
C, 


D 


P  (N) 

Cj.  - 


'M 


'N 


'"SR 

£ 

^SR 

g 

H 

h 

Ir.  Ip 


K 

k 

L 

^TR 

n 

P 

Q. 


Frontal  area 
Drag  coefficient 

Newtonian  drag  coefficient 

Computational  constants 

Coefficie^it  of  moment  about  pitch  axis 

Coefficient  cf  moment  about  yaw  axis 

Specific  heat 

Stiffener  ring  depth 

Modulus  of  elasticity 
Safety  factor 

32.284  ft/sec^ 

Total  enthalpy 

Specific  enthalpy 

Roil,  pitch  moment  of  inertia 

Pitch  moment  of  inertia  with  respect  to  plane  normal  to 
the  X  axis  through  the  origin. 

Experimental  constant 

Computational  constants 

Coefficient  of  thermal  conductivity 
Latent  heat  of  sublimation 
Transition  parameter 

Number  of  computational  intervals  in  a  given  section 
Pressure 

Total  heat  transfered  to  the  body  at  any  point 
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Average  cold  wall  heat  absorption  capability 
Heat  flux 

Heat  flux  at  tinne  t 

o 

Radius 

Reynolds  number  based  «,  i  laminar  boundary  layer  momen¬ 
tum  thickness 

Radial  coordinate 

Temperature 

Maximum  permissible  temps,  ature  on  inside  of  wall  at  end 
of  ablation 

Initial  temperature 
Ablation  temperature 
Time 

Total  heating  time 
Velocity 

Norro.alized  velocity  constant 

Dimensionless  ablation  rate 
Weight 

Longitudinal  body  position  coordinate 
X  coordinate  of  center  of  gravity 
X  coordinate  of  center  of  pressure 

Total  wall  thickness 

Altitude 

Coordinate  normal  to  body 
Thermal  diffusivity 

Scale  height  in  atmospheric  density-altitude  relationship  = 
4.5  X  10-^/ft 

Ratio  of  specific  heats 

Ballistic  coefficient  =  W/Cj^A 
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Boundary  layer  thickness 
Body  angles 

Re-entry  angle  of  trajectory  with  horizontal 

Coefficient  of  viscosity 
Density 

Reference  density  in  atmospheric  density -altitude 
relationship 

Thickness /depth  ratio  of  stiffener  ring 

Material  thickness 
Body  angle 
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Subscripts 

B 

E 

L 

M 

N 

o 

P 

P/L 

q 

t 

w 

a 

5 

CO 

CO 
C.  W. 

H.S. 

«-S-abl. 

®‘RES. 

INS. 

P/lu 

SR 

ST 

TOT 

TR 

R 


Base 

Re-entry 

Laminar 

Material 

Nose 

Stagnation 
Cover  plate 
Payload 

Partial  derivative  with  respect  to  roll  velocity 

Turbulent 

Wall 

Partial  derivative  with  respect  to  angle  of  attack 

Edge  of  boundary  layer 

Free  stream 

Cutoff 

Cold  wall 

Heat  sink 

Heat  sink  ablated 

Total  heat  sink 

Heat  sink  residual 

Insulation 
Payload 
Stiffener  ring 
Structure 
Total 

Transition 

Reference 
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I.  INTRODUCTION 

The  purpose  of  this  report  is  to  serve  as  an  aid  in  running  the  IBM  704 
Re-entry  Vehicle  Synthesis  Program  for  parametric  design  studies  of  re¬ 
entry  vehicles.  The  report  is  essentially  divided  into  two  parts:  the  first 
part  presents  the  basic  equations  that  are  utilized  and  the  assumptions  that 
have  been  made,  while  the  second  (sections  IX  and  X)  deals  with  the  actual 
computational  procedure  that  the  machine  uses  to  solve  the  equations  for 
the  first  part  in  the  case  of  a  given  re-entry  vehicle.  A  straight  line  trajec¬ 
tory  is  assumed  and  the  heating  «inalysis  is  for  the  case  of  a  non- reacting 
cold  wall,  with  heat  conduction  into  the  wall  and  boundary  layer  transition. 
The  computer  output  consists  of  values  of  ablation  and  insulation  thicknesses, 
vehicle  weight,  and  aerodynamic  coefficients.  All  of  the  heat  equations  were 
integrated  over  the  trajectory  prior  to  programming  so  that  the  computer 
output  consists  only  of  total  rather  than  step-wise  values. 
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II.  TR.AJECTORY  EQUATIONS^ 

The  following  trajectory  equations  were  derived  by  Allen  and  Eggers^ 
under  the  assumptions  of  constant  drag  coefficient,  negligible  gravitational 
effects,  and  an  exponential  variation  of  atmospheric  density  with  altitude: 


oo 


=  u_  exp 


,  -PY  1 

2P 

W  ^ 

sin  0^ 
E 

where 


-  -(3Y 
p  =  pe  ^ 


du 


00 


dt 


2  E 

-  p  sin  0  u^  In - 

E  oo  u. 


CO 


oo 


e  ^  “oo 


III.  THERMODYNAMIC  ASSUMPTIONS 


(1) 


(1) 


(2) 

(3) 


a)  The  air  passing  over  the  vehicle  is  a  dissociated  gas  in  thermal 
equilibrium . 

b)  The  conditions  at  the  edge  of  the  boundary  layer  at  any  point 
along  the  body  can  be  computed  from  the  stagnation  point  con¬ 
ditions  by  means  of  an  isentropic  expansion,  since  the  energy 
losses  to  the  body  are  relatively  small. 

c)  The  viscosity  varies  as  the  square  root  of  the  enthalpy. 

d)  An  average  value  of  y  =  1.2  will  be  used. 

Under  these  assumptions  the  following  relationships  hold: 
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IV .  LAMINAR  CONVECTIVE  HEATING^  ^  ^ 
The  basic  equation  to  be  solved  is 


q 


L 


k 

w 


(4) 


(5) 


(6) 


(7) 


This  equation  can  be  transformed  by  using  the  Levy  transformation 
combined  with  the  Mangier  transformation  for  bodies  of  revolution.  The 
following  assumptions  are  made: 

a)  The  variation  of  heat  transfer  rate  with  local  pressure  gr.adients 
is  negligible. 

b)  pp  may  be  assumed  constant  across  the  boundary  layer  if 
evaluated  at  some  average  enthalpy. 

c)  The  free  stream  static  enthalpy  is  negligible. 

Under  these  assumptions  the  laminar  cold  wall  convective  heating 
rate  can  be  written  as 


C,  u 

1  00 


1/2  .  j  _  1/2 


■B 


(8) 


ST  L/TR -60- 0000 -09 103 

Page  4 


where  Fj(r/Rg)  represents  the  variation  of  local  heating  rate  .^round 
the  vehicle  and  is  a  function  of  geometry  and  pressure: 


The  limit  of  this  expression  at  the  stagnation  point  r  =  0 


(9) 


is  found  to 


be: 


F^{0) 


1/2 


(10) 


By  equation  (6)  can  be  expressed  as  a  function  of  u^.  An 
approximate  expression  for  Pq/Pqq  determined  as  a  function  of  ^ 

from  gas  tables: 


c  oo 


(li) 


From  the  above  relationships  and  equation  (2),  the  time -dependent 
variables  in  equation  (8)  can  be  expressed  in  terms  of  u^-  Integration  of 
the  resulting  equation  from  the  initial  value  of  u^  =  u^,  to  any  arbitrary 
value  of  u^  gives  the  total  laminar  convective  heat  transfer  per  unit 
area  as 


lli 


1/4 


^B/ 


1/4 


1 


R 


B 


u. 


erf  '\/2.  1  52  In 


u 


(12) 


oo 


where 

ir  _  i/"' 

■  ‘L  ”  "E 

and  is  an  experimentally  determined  constant 


^  2.  152  ,  .  „  .-1/2 

u,-,  VSin  »_) 
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V.  TURBULENT  CONVECTIVE  HEATING^ 

In  this  case  the  basic  equation  to  be  solved  is  the  integral  form  of 
the  energy  equation: 


d 


Is  r  rpu 

(  Jo 


(H  -  Hg)  d 


-  r 


(13) 


The  wall  heating  rate  may  be  computed  from  the  Blasius  turbulent 
flat  plate  solutions  in  the  form: 


H„) 


(14) 


where  Re^  is  defined  as 


Re. 


pu6 


(15) 


It  is  assumed  that  pu  and  H  vary  with  y/6  as  follows: 


JPiL 


H  -  H 


H  ,  .1/7 

rf-ii) 


(16) 


This  yields  an  expression  for  6  which  can  then  be  substituted  into 
equation  (13)  to  give  the  turbulent  convective  heating  rate.  If  the  starting 
length  is  measured  from  the  nose,  the  turbulent  heating  rate  can  be  written 
as: 


^t  = 


C3p'u(p')^^^(H5-  H^) 


5/4  ,  ,  ,  J/4  , 

p  u  (p  )  dx 


TT73- 


(17) 


t  ( 


If  P  .  p  are  the  values  of  p,  p  evaluated  at  some  average  enthalpy 
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to  include  compressibility  effects  across  the  boundary  layer,  and  the  ratio 
1 /4  ;  1 /4 

p'  (p.’)  /pp  IS  averaged  along  the  body,  the  equation  {i7)  becomes: 


u 


y  F 


(18) 


(19) 


As  in  the  case  of  laminar  heating,  the  quantities  p^,  h^,  p  ^  can  be 
expressed  in  terms  of  u^,  and  the  resulting  equation  for  q^  integrated 
to  yield  an  expression  for  the  total  turbulent  heating  per  unit  area; 


=  K. 


1  F 

/  r  1  ■ 

1^)  ^ 

HI 

r  (4/5)-r  (4/5) 

A 


(20) 


where 


K 


u, 


2.  485 


(sin  0g) 


1/5 


and  K'^  is  an  experimentally  determined  constant  (4/5)  represents 

the  incomplete  gamma  function  with  X  =  2.  4853  In  Uj,/u 


It  should 


TR 


be  noted  that  there  is  a  stronger  dependence  on  (W/C^A)  in  turbulent  flow 
than  m  laminar  flow. 
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VI.  TRANSITION  FROM  LAMINAR  TO  TURBULENT  FLOW 

The  transition  criterion  selected  for  this  program  is  based  on  the 
momentum  thiclcness  Reynold's  number  which  is  given  by 


Substituting  (4),  (5),  and  (6)  into  (21)  gives 


(22)  may  be  re-written  as 


(24) 


Assuming 


0.  36 


(25) 

(26) 


I 
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(24)  becomes 


/R, 


0.  36 


'6 


1^3/ 


=  C 


7  a. 


P  ^ 
^oo  oo 


0.  88 


(27) 


Substituting  (3)  into  (27)  gives 
2 


R  ' 
®0 


\  ^3/ 


Cg  sinO^ 


W  1 

1  u„) 

fv\ 

1  “j 

0.  88 


u. 


0.  6  0.  28 


u. 


(28) 


Defining 


(28)  becomes 


N  = 


“e' 

00 

\  I 


/U 


.0.  88 


on 


\  / 


(29) 


N  = 


/R  ' 


3i 


w  1 

1 

c  1 

1  .  “e  1 

sin 

»e1 

l^D^i 

9  I 

[23,  477] 

0. 

35/ 

18.  2 


(30) 


or,  at  transition, 


TR 


N 


TR 


ir 


3  / 


'N 


TR 


(31) 
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where 


'N 


=  C, 


TR 


u 


0.  28 


E 


123,  477 


sin  1 

TTlfj 


<Po^Poo>j 
18.  2 


(32) 


The  reference  enthalpy  and  viscosity  are  given  by 


ho  =  2,  340,  000  ft-lb/slug 

^  6 

=  0.  551  X  10  slug/ft-sec 


Then 


Cg  =  60,  000  ft/lb 

Only  one  branch  of  (29)  is  entered  in  the  program: 

(u  /u„)  =  1.  At  (u  /u„)  =  0.  32,  a  peak  of  N  =  0.  41802  occurs,  hence, 
00  E  00  E 


there  is  no  u 


for  body  stations  at  which  Nrpo  ^  0.  41802,  indicating 
TR 


that  only  laminar  heat  transfer  is  experienced  at  these  stations.  To  prevent 
the  possibility  of  aft  stations  becoming  laminar  (through  the  mechanical  use 
of  (31)  once  a  more  forward  station  has  become  turbulent,  is  com¬ 


pared  with  the  value  of  at  the  preceding  station  (N^^). 


If 


N 


TR 


=  then  the  value  of  is  used  to  force  turbulence  at  least  as 


'TR 


soon  as  it  occurs  at  the  immediate  forward  station.  If  N 
;he  program  uses 

VII.  HEAT  SHIELD  AND  INSULATION  THICKNESSES^ 

The  basic  Fourier  heat  conduction  equation  iS  used: 


TR 


N 


TR’ 


then 


2 

a'^T 


1  aT 
Q  TT 


(33) 


and  it  is  assumed  that  the  heat  shield  is  a  semi-infinite  slab.  Ablation  is 
accounted  for  by  the  designation  of  a  moving  surface  in  the  boundary  condi¬ 
tions.  It  can  be  shown  that  for  a  heat  input  to  the  surface  of  the  form 
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q  =  q 


m 


>1" 


(3^) 


the  following  equation  holds  for  the  normalized  velocity  of  ablation 

il/2 

I  w  ' 

V"  =  V  ‘  " 


(35) 


where  V  is  determined  from 


yT(i  -  v^)yH-rq^ 

MT3  -1-^1 


exp 


/.2 

,r2 

/  q  1 

V 

^  0 

0 

2 

7T 

\^Pm 

L 

erfc 


\ 


q  V 

^m  o 


(36) 


The  quantity  ^ound  by  integrating  equation  (34)  over  the 

total  time  of  heating  to  yield 

(vv^:) 

where  is  determined  from  the  procedure  of  the  preceding  section  and 

t^  is  the  total  heating  time.  Once  known  it  is 

then  possible  to  solve  for  the  thickness  of  material  ablated,  Ttt  c 


If  the  heat  shield  and  insulation  have  similar  thermal  properties,  the 
total  wall  thickness  prior  to  ablation,  x.^,,  can  be  found  from  the  following 
equation: 


STL/TR-60-0000-09103 
Page  1 1 


T.  -  T 
1  o 

T'"-T' 

s  o 


erfc 

Xt  ^ 

J 

erfc 

j 

(39) 


If  the  thermal  properties  of  the  heat  shield  and  insulation  are  different 


and  if  the  insulation  cannot  exceed  some  given  temperature,  T. 
the  total  heat  shield  thickness,  t 


then 


H.S. 


can  be  determined  from 


max 


T.  -  T 
1  o 

max 

- 

s  o 


^H.S.  \ 

erfc 

maxj 

1  / 

^H.S  \ 

erfc 

ABLJ 

j 

(40) 


The  total  insulation  thickness, 


INS.  ’ 


is  then  given  by 


■^INS. 


H.S. 


max 


f^INS.  \ 

) 

rH.S.) 

1/2 


(41) 


where  it  is  assumed  that  the  p.,  C-.,  product  of  the  two  materials  is 
constant. 

Vm.  VEHICLE  GEOMETRY 

The  vehicle  shape  is  defined  as  follows  (see  Figure  1); 

a)  A  spherical  segment,  either  solid  or  hollow. 

b)  An  abritrary  number  of  (at  least  one  but  not  more  than  ten) 
conical  frustums. 

c)  A  cover  plate  at  the  rear  of  the  vehicle  which  may  be  either 
convex  (dashed  line)  or  concave  (solid  line) 
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The  number  of  computational  intervals  in  each  section  of  (a)  and  (b) 
can  be  arbitrarily  chosen  from  one  to  eighteen  The  number  of  frustums 
and  the  following  quantities  must  be  specified  to  determine  the  geometry 
of  the  vehicle. 

must  always  be  sp'  't  ed.  If  0  j  is  not  given  it  will  be 
assumed  equal  to  (tt/Z  -  which  is  the  tangent  case. 

t>)  is  always  specified  but  in  any  other  (i^^)  section  any  two  of 

the  quantities  [0  ,  Ax.,  R..,)  maybe  selected,  x.  ,  and  R. , 

\  t  lid)  1 1  1 1 

are  computed  by  the  program. 

c)  ^2  is  assumed  specified  and  Rp  =  f  (^2'  ^3^’  Provision  is 
made  for  the  selection  of  Rp  and  the  determination  of 

^2  ~  ^  (Rpi  Rg)-  When  either  or  Rp  are  given  as  nega¬ 
tive  the  cover  plate  will  be  defined  as  convex. 

If  the  nose  cone  is  specified  as  solid  it  will  be  assumed  to  be  heat 
sink  material,  with  the  radius  of  the  spherical  segment  equal  to 

S.  ^INS.  ^  juncture  of  the  first  section  plus  t  if  the 

structure  is  on  the  outside. 

All  the  point  radii  will  be  normalized  by  R^.  it  will  be  assumed 
that  the  vehicle  geometry  does  not  change  for  consecutive  runs  unless 
specified  and  then  only  as  specified. 

IX.  FIXED  WEIGHT  AND  PAYLOAD 

An  arbitrary  amount  of  fixed  weight  and  payload  can  be  defined  in 
each  section  The  center  of  gravity  (x)  of  the  payload  is  measured  in 
inches  from  the  front  of  the  frustum  in  which  it  is  located;  the  x  of  the 
fixed  weight  is  taken  to  be  at  the  x  of  the  section,  independent  of  the 
distribution  of  the  fixed  weight.  The  pitch  and  roll  moments  of  inertia 
of  the  payload  are  given  about  the  payload  center  of  gravity;  the  Lp  and 

Ir  of  the  fixed  weight  are  given  about  the  section  center  of  gravity  deter¬ 
mined  without  the  fixed  weight. 
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Axial  symmetry  is  assumed  for  the  entire  vehicle.  It  will  be  assumed 
that  these  values  do  not  change  for  consecutive  runs  unless  otherwise 
specified. 

X.  MACHINE  COMPUTATIONS 

This  section  and  Section  Xl,  describe  the  computation  procedure  of  the 
machine.  The  various  sections  of  the  total  computation  are  listed  below. 

These  should  be  used  with  the  corresponding  flow  diagrams  given  in  Section  XI 
to  understand  the  procedure  of  any  given  section  of  the  computation.  The 
sections  below  are  listed  in  the  order  in  which  they  occur  during  the  compu¬ 
tation  and  the  chain  of  computations  is  shown  on  the  first  diagram  in  Section  XL 

A.  SHAPE  (Body  geometry  calculation) 

This  portion  of  the  program  determines  the  geometry  of  the 
vehicle  from  the  inputs  described  in  Section  VIII.  The  coordinates  of  the 
points  (x/Rg,  r/Rg)  are  computed  and  stored  in  the  point  bins.  At  present 
the  program  divides  each  section  into  ten  intervals.  No  flow  diagram  is 
included  since  the  computation  is  quite  straightforward  because  of  the  simple 
geometries  involved.  A  complete  set  of  applicable  flow  diagrams  is  con¬ 
tained  in  Section  X. 

B.  PRESS  (Pressure  distribution  calculation) 

The  pressures  on  the  vehicle  are  supplied  to  the  program  by 
means  of  linear  interpolation  from  tables  of  pressures,  for  Mach  20  and 
altitudes  of  100,  000  ft  and  65,  000  ft.  For  the  first  two  sections  these  are 
pressure  values  for  sphere-cone  combinations  that  have  been  determined 
by  a  method  of  characteristics  solution,  and  are  an  integral  part  of  the 
program.  For  the  remaining  sections  the  pressure  tables  are  a  part  of 
the  input.  If  desired,  a  Newtonian  table  of  pressures  which  is  provided 
in  the  program  may  be  used  for  the  remaining  sections. 

C.  FFF  (Body  function  calculations) 

The  expressions  to  be  computed  are  F^,  F^,  F^,  as  given  by 
equations  (9),  (19),  and  (23).  The  computational  procedure  is  outlined  on 
the  included  flow  diagram. 
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D.  CPM  (Cover  plate  calculations) 

The  following  quantities  as^sociated  with  the  cover  plate  are 


computed; 


X  =  X 


B 


Rp  ^ 

-y-  (1  -  cos  =  - 2 - 


W 


2it 


2ir 


ST 


3  ^ST  ^ST  ^P  "  73  ^ST  "^ST  ^P 


12 

W  „2 


12' 


^B  '  "‘C. 


3^  Rp  (2  -  cos  0  -  cos  0)  .  — 2 


I>2 

X  ~  3g  772  ^  — T" 


12 


12  g 


x^  -  (x^  +  Rp  -  x)^l 


All  densities  are  given  in  Ib/ft  .  Lengths  are  given  in  inches  and 

weights  are  given  in  pounds.  The  moments  of  inertia  will  have  the  units  of 
Z 

slug-ft  •  The  computational  procedure  is  outlined  on  the  included  flow  diagram. 

E.  ARDC  (Aerodynamic  coefficient  calculations) 

The  following  Newtonian  aerodynamic  coefficients  associated  with 
the  vehicle  are  computed; 


'D 


'R 


3.  92 
^B 


B 


0 


1^1 
o  , 


'D 


(N)  _ 


3.  92 


B 


sin^'  0  dr 


'N 


R 


3  92 

Jo 


B 


R 


B/ 


cos  0  dr 
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3  92  R 


B 


"•^B  1 

f  ^ 

- 

1  ) 

- 

r 

i 

! 

sm  0  +  j 

. 

1  1 

cos  0 

Jo 

(cos  0)  dr 


C  P. 


M 

a 

7Z - 


N 


B 


a 


Since  these  are  straightforward  calculations  no  flow  diagram  is 

included 

F.  ITC  (Iterative  stagnation  pressure  calculation) 

The  equation  currently  being  used  for  this  program  is 


P 

o 


3.8  X  10 


-9 


W 


CdA 


u^  sin  0, 


where  u^  =  23,477  ft/sec  and  0^  =  20.5°.  This  value  of  P^  corresponds 
to  P^  at  maximum  dynamic  pressure . 

G.  STRG^"^^  (Structure  calculations) 


The  structure  thickness  the  structural  weight 

stiffener  ring  depth  and  stiffener  ring  weight  are 

This  section  of  the  computation  considers  only  the  hollow  nose; 
of  the  solid  nose  is  included  later. 


«ST' 

computed . 
the  case 


For  the  nose: 


“^ST  =  ^ 


P 

o 


1/2 


W 


2ir 


ST 


12 


PcT  ^ 


3  ^ST  ST  "N  12 
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SR 


2!sR  ^  j,3 

E  12  12 


.1/4 


W 


2ir 


SR 


SR 


12 


_  j 

3  "’SR  ^ST  SR 


R 


12 


SRi 

T 


For  the  frustums: 


ST 


W 


1  121 
cos  "5 

TT 


(Pi'  / 

1^:)  1 


X. 


i2 


^il 


R.  .  +  R.  , 
1 1  iZ 


ST 


12  cos  0 


PST  V  (^2  -  =‘ilj  (^1  "  ^-iz) 


SR 


f 


SR  P 


'SR 


1/4 


W 


SR 


2ir  ,2  1„  -“SR 

^SR  PsT  ^SR  ri2  ■  ~2 


Each  stiffener  ring  joins  two  consecutive  sections,  including  one 
joining  the  nose  to  the  first  frustum  and  one  joining  the  last  frustum  to  the 
cover  plate . 

In  the  above  equations,  f  is  a  safety  factor  which  is  assumed 

Siv 

to  be  1.25  for  all  computations.  The  modulus  of  elasticity  E  may  be 
varied  from  section  to  section.  The  is  a  thickness /depth  ratio, 

^SR^^SR’  stiffener  rings.  The  pressure  P  is  taken  to  be  the 

maximum  value  of  pressure  in  the  section  under  consideration. 
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H.  THSN  (Cold  wall  heat  sink  thickness  calculation) 


In  this  section  the  value  of  t„  ...  /R„  is  computed  as 

O  .  W  Jl3 


^  ^H.S. 


This  implies  only  cold  wall  heating  with  no  conduction  of  heat 
into  the  body.  The  ^  thus  found  differs  from  the  value  of  Tj^  g 
given  in  the  output.  By  a  later  sophistication  of  the  program  only  the 
value  of  Q.p  (total  heat  transferred  to  the  body  at  any  point)  is  computed. 
This  value  is  then  used  as  input  to  equation  (37), 


In  the  computation  of  the  criterion  used  to  determine 

whether  heating  is  laminar  or  turbulent  is  the  value  of  the  transition 
parameter  defined  as 


N. 


TR 


1  N 


TR 


R, 


'0 


TR 


iTHTr^ 


®TR 


is  given  as  input  for  each  case  and 


is  defined  as. 


6,  00  X  lO'^ 


“e 

0.28 

[Sin  e^i 

1^0  ^ 
i‘’“i 

E 

23.477] 

l~OF~i 

18.2 

where 


,/p 


IS  a  tabulated  linear  function  of  u^  Thv'  equation  for 


(^o  ^  ^ 

Q.p  (see  flow  diagram  and  N .  T  C.  calculation  sheet)  is 
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Qt  =  Q  Ph  s  (H.S.B.)  d  /(HSC)  erf  (e)  +  e 


1  41 

|4i' 

r  k  r 

xco' 

r  -r 

xtr'  ^  [ 

where 


i  =  2.152  In 


u. 


1/2 


CO 


This  equation  is  modified  for  the  various  values  of  as  follows: 

1  ix 

for  N  ^  0.4179,  set  erf  (^)  =  land  r  (4)  =  F  (41 

Xtr  Xr.nw/ 


^TR  “  ^TR  ^TR 


^TR  ^TR  ^TR 


There  are  tables  stored  in  the  computer  giving; 

u  1 


-  f  (N^j, 

erf  (4) 

=  f 

/U 

oo 

^  v 

Xtr  >5; 

=  f 

u 

oo 

I.  TIN  SR  (Heat  sink  and  insulation  thickness  calculation) 
For  this  computation  equation  (37)  is  interpreted  as: 

♦ 


Q 
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from  which  is  found.  The  value  of  in  this  computation  is  currently 

given  in  the  program  as  30  seconds.  This  value  of  substituted  into 

equation  (36)  giving 

....  .2 


V“H.  S.  ’’cwPh.S. 

2A‘'  'T=  -  V 


exp 


Q  T 

_  cw 


or,  in  terms  of  the  con.stants  defined  on  the  flow  diagram, 


K3  (1  -  V^)  exp  (K^  V^)  erfc  =  1 


V 


With  the  introduction  of  the  new  variable  V  =  V  K, ,  this 

o  o  1 


equation  becomes; 


*  K-  .  .2 

g  (V^)  =  (K,  -  V^)  exp  V 


1 


T 


erfc  V  -1=0 
o 


which  is  the  equation  utilized  by  the  computer.  If  V  =*3  the  quantity 

jjc 

erfc  must  be  computed  by  the  complementary  error  function  routine 

on  the  computer  for  accuracy.  If  V  •<^3  the  faster  computer  route  for 
❖ 

erf  may  be  used  and  the  result  subtracted  from  1. 

*  >(' 

The  value  of  V  for  which  e  (V  )  =  0  must  now  be  found.  It 
o  ^  '  o' 

is  easily  seen  that  g  (0)  =  -  1  and  g  (Kj^)  =  -1,  so  that  the  required 

value  of  V  lies  between  0  a.nd  K,  if  K,  =•  1.  A  binary  search  is  con- 

o  ^  ^  ^ 

dlucted  for  the  root  of  g  (V  );  the  V  selected  is  that  for  which  the 

absolute  value  of  the  difference  between  consecutive  values  of  V  is  less 

o 

than  0.  01. 


The  thickness  of  heat  sink  ablated  is  then  computed  from 
equation  (38).  Then  equation  (39)  can  be  solved.  If  77  8.  5  the  value  of 

^INS  extremely  small— as  indicated  by  equations  (39),  (40),  and  (41)— 
and.  '''jjvjg  is  set  equal  to  the  minimum  specified  value  of  0.  1  inch.  The 
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computational  method  to  determine  erfc  r|  is  analogous  to  that  described 
above  for  erfc  .  if  the  ablation  temperature  is  equal  to  or  less  than  the 
maximum  permissible  inside  wall  temperature,  the  computed  value  of 

^INS  zero;  thus  is  again  set  equal  to  the  minimum 

specified  value  of  0.1  inch. 


In  the  Section  of  the  flow  diagram  enclosed  in  the  dashed  lines  a 
binary  search  is  conducted  for  which  is  the  total  wall  thickness  for 

homogeneous  material  construction,  provided  the  computed  >  0.1  inch. 

This  quantity  is  found  just  before  the  flag  check  is  entered  for  the  first  time. 


The  program  then  proceeds  out  of  the  flag  check  on  the  path 
marked  (1)  and  checks  to  determine  whether  T:  =  T  If  qn  the 

insulation  and  heat  sink  are  of  the  same  material  and  t  is  computed 

as  shown.  If  not,  the  heat  sink  and  insulation  materials  are  different. 

After  the  flag  is  changed  to  show  that  x,p  routine  is  completed  and  the 
changes  noted  on  the  flow  diagram  are  made,  the  procedure  enclosed  in 
the  dotted  lines  is  again  utilized  to  compute  the  total  heat  sink  thickness, 

■^H.S-n.ax'  ^INS.  computed  as  shown  for  the  case  of  different 

insulation  and  heat  sink  materials. 


J (Heat  sink  and  insulation  weight  calculations) 

In  this  section  of  the  computations  the  weights  of  the  heat  sink 
and  insulation  materials  for  a  hollow  nose,  solid  nose,  and  the  frustums 
are  determined. 


For  the  solid  ncse: 


Ztt  p 

'  ‘'OS 


W 


INS. 


0 
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where  t  =  +  x  „  at  the  juncture  with  the  first  section  plus  x  „ 

IN  S  H .  o  •  ^ 

if  the  stricture  is  on  the  outside  of  that  section 

In  the  case  of  the  hollow  sections  the  following  computational 
technique  is  used- 


'ms.  =  to  compute- 


’'ins.  *  '"h.s 


=  T.,  to  compute  — 
K  P, 


Then  the  heat  sink  and  insulation  weights  are  determined  from 
the  following  relationships. 


'"ins  =  ■  Pins.  '  77 


'"2  '"1 


«  K  s  ■  T 


For  the  hollow  nose: 


3  sin  <pd  p 
3 


where  dp  =  p^/n  in  radians,  and  is  replaced  by  (R^^.  +  x^.^)  if  the 

structure  is  on  the  outside  of  the  nose. 


For  the  frustums: 


W  ^2 

Pxx  "  ""  B 


cos  6 


r  dx  =  2 


2  r^r  ""k  ' 

J  +  2  cos  ej 
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where 


ds 


^2  -  ^ 
n  (cos  6) 


K.  WTOTS  and  WTOTB  (Iterative  ballistic  coefficient  calculation) 


These  two  operations  are  straightforward  and  require  no  flov; 
diagram.  The  computation  equations  are: 


"tot  = 


UNO. 


where  j  denotes  a  section  of  the  vehicle; 


Atot 


i  B  cos  0  cos  6 


This  expression  will  also  include  a  “Tgrp  term  if  the  structure 
is  on  the  outside. 

The  ballistic  coefficient  is  then  computed  as 


^i  +  1 


W 

TOT 
*  ^TOT 


L.  TEST 

This  section  of  the  computations  is  explained  on  the  flow  diagram 

M.  XIXIR  and  IRXIX  (Moment  calculations) 

The  equations  used  for  the  computations  in  these  sections  are 
given  on  the  flow  diagram.  The  primed  quantities  refer  to  volumes  rather 
than  masses;  that  is,  I  refers  to  the  volume  moment  of  inertia  and  M* 
to  the  volume  first  moment.  The  XIXIR  computations  are  only  for  the  nose 
and  are  followed  directly  by  the  IRXIX  computations  for  the  frustums. 
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N.  MSTR 

This  set  of  computations  is  adequately  described  by  the  flow  diagram 

O.  TBM  (Tabuiati'^n  computations) 

These  computations  are  straightforward  and  a  flow  diagram  was 
not  considered  necessary. 


TOT 


2  w.  (X). 


TOT 


TOT 


= 

j 


11  W.  r 

1 

I  +  -£  + 

+ 

1  tot  j)  j 

^  FIXED  ^P/L 


W 


P/L 


12^  g 


^TOT  '  ^P/Ll. 


W 

2  ^  FIXED 


12  g 


XpoT  ' 

J 


‘R 


TOT 


weight . 


'  2  (^1- 
J  ^  'j 

These  quantities  are  computed  both  with  and  without  heat  sink 


P.  ETC  (Miscellaneous  calculations) 

The  following  Newtonian  stability  parameters  are  also  computed 
by  the  program  at  the  end  of  the  other  computations: 

nR 


3.92  •  R 


B 


'M 


B  Uq 


B 


(r/Rjg)  cos  2  6  ((r/R^)  tan  8  + 


B' 


;c 


l^B/ 


TOT' 
^B  j 


>dr 


Stab,  Par 


'N, 


D 


C  ^  C 
m  m 

q  a 


TOT 

B 

L-X  . 

12 

TOT 


^  ^  ^B 
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XI.  FLOW  DIAGRAMS 

This  section  contains  the  following  flow  diagramst 

1.  General  Control  for  Computations 

2.  FFF  Calculations 

2a.  FFF  Calculation  Notes 

3.  CPh'  Calculations 

4.  N.T.C.  Calculations 

5.  ST RG  Calculations 

6.  THSN  Calculations 

7.  TINSR  Calculations 

8.  Test 

9.  XIXIR  Calculations 
10  I RXIX  Calculations 
11.  MST R  Calculations 
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SETUP  -  SCUP" 


Diagram  1.  General  Control  for  Computations. 
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Diagram  2.  FFF  Calculations. 
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STARRED  QUANTITIES  REPRESENT  VALUES  AT  THE  POINT 

UNDER  CONSIDERATION  -  UNSTARRED  QUANTITIES  ARE  FROM 

PRECEDING  POINT. 


G  (!) 


'■/Rb  21±i  r 


o 


ds 


r/RB 

G(2)  =y 


y+7 

8y 


Hi 
y 


ds 


FRUSTUMS 


SINO 


Diagram  2a.  FFF  Calculation  Notes. 
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4  /  “e  \0-28sfNeE  E 

6.00X10  (^23477j  0.35  ’  18.2 


WHERE  {Po/P<d)  is  tabulated  VS.  IF  hg  =  250,000 
—  OTHERWISE  GO  TO  ERROR  AND  STOP. 


(2)  HSA  :  0.748 


f-jLJ 

\23477 ) 


0.335  /, 


Ro\  0-3 


'SIN  ^ 

0.35  ■  12  j 


(3)  HSB=  (R^/Rb) 

/  r  v/4 

(4)  HSC  - 


(5)  CQL  I  899 


/Ul. 


,2.15  , 

''E  \  fs\u4>\ 

23477/  \  0.35  ) 


Diagram  4.  N.  T.  C.  Calculations. 
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NOTE 

/7^.  A  »  J 


P  RT 

VALUES  G"  ^cT  oStu  ARt  equal  TO  — ^ 

12^ 


Diagram  5.  STRG  Calculations. 


Diagram  6.  THSN  Calculations 
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note: 

i'  :  i-t-l  IN  PRECEDING  EQUATION 
A  =  005  IN  CURRENT  PROGRAM 
A*=  0.01  IN  CURRENT  PROGRAM 

THE  EQUATION  BEING  USED  F0R"ITC"  IN  THE  CURRENT  PROGRAM 
is;  Pq  =  3.8  X  io'^Au^e  sin  % 

—  THE  MACHINE  THEN  PROCEEDS 

TO  THE  "STRG"  CALCULATIONS 


Diagram  8.  Test. 


FRAMES 


fc»— - - -  ■ 


SET  COMPUTER  TO  NOSE 


rsTCO  ? 


c 
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2rrpHs  (^X  " S  A8L  )' 1''°^  '  0  ({\  +  COS  <p) 


'2  '^INS  +HSRES 


V 


r 

(_  “  ^HS  ABL  J  3 


Pus  (COS  4>-\ ) 
'2^  W^^S  ABL 


„  sfc  +  cos  4>)  „ 
Rx  - g -  Rx 


3 

-(Rx  Ths  aBL  '> 


(!  +  COS<^) 
8 


(Rx  ^HS  ABL 


.  R  ,  cos  (rx  -  -l^)]  -  h  -»  s  abJ’ [^  -( " X  -  ' « s  abl)  cos  ([ ".? 
'msabl)  (  T  +(''x-'-HSAeL)cos+((-^^  i^-^C05#  -  Jj| 


Rni 
■  2 


—  +  2 


5  r_  >l5l 

X  "  *■  hS  ABL;  I 


R  X  “  R  V  “  r 


s^4> 


+  2  <  Rx  ■  HS  ABL 


_L 


COMPLE'^ES  NOSE  -  GO  '^O  ''RST  fRUSTuV  'RXIX 


Diagram  9.  XIXIR  Calculations 
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Diagram  10,  IRXIX  Calculations. 
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note; 

0  throughout  calculation  use  (IA4  9)  for  g 
@  although  the  moment  values  for  the  rings 

IN  THE  SOLID  NOSE  ARE  COMPUTED, THE Y  ARE  NOT  uSED 


Diagram  11. 


MSTR  Calculations 
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